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ABSTRACT 


Current methods of shock design and analysis rely primarily 
on the input of lerge amounts of experimental data, assembled from 
underwater explosion tests, to describe the shock motion inputs to 
their systems. The object of this investigation is to arrive at a 
more complete characterization of the shock motion input, a net 
effect representing the interaction between the water applied 
pressure loading, the hill structural response and the equipment - 
foundation response. 


The model chosen is the cross section of a simple 
cylindrical shell with an inboard mounted simple oscillator, 
immersed in an acoustic medium. The input to the water-shell 
system is an exponentially decaying plane step wave incident 
normal to the cylinder longitudinal axis, 


The analysis employs the use of four-pole mechanical 
impedance techniques and Laplace transform methods. An analytic 
expression is derived, which has the form of a mechanical cuad- 
ripole impedance. It relates shock wave pressure to equipment 
response, as a function of sheli geometry aad oscillator mass and 
natural frequency. The lack of rigor associated with the description 
of a Gistributed mass system as a point phenomenon along with other 
restrictive conditions, cause the expression to be better termed 4 
"“quagi-impedance” . 


The problem formulation end solution has limitations on 
poth the mathematical and physical side and at this stage the re- 
sults are to be regarded mainly as exploratory. The " quasi-inpedance" 
4s now capable of predicting qualitatively the effects of parameter 
variation on the response of simple equipments to shock wave inputs. 
Additional refinements of the shell model and comparison with experimental 
data will enable a more rigorous evaluation of quantitative effects. 
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ratio of the longitudinal length of shell 
feeling the influence cf the equipment reaction 
to the actual equipment length 

mass of shell per unis area 

modal ox generalized masses of shell per unit 
area 


mass of equipment model 


mode amber and equivalentiy the number of 
circumferential waves of mode 


free-field shock wave pressure 

radiated and diffracted wave fluid pressure 
generalized or modal components of P,, P. 
shock wave peak pressure 


generalized or modal deflection in the radial 
direction and its derivative with respect to 
time, the generalized velocity 


Laplace transform of a 


cylindrical coordinates 
complex frequency operator 


roots of characteristic equation 


tine 

non-dimensional time = ect/a 

step function 

free-field shock wave particle velocity 
radiated and diffracted wave fluid particle 


velocity 
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= modal components of Uys U,. 


V i = point velocity where 1 refers to the particular 


component output or input point 


". = model component of V, 
Yo = peak velocity of qasi-~impedance functioa 
J = vadlelly inward net velocity of shell as a 


X(s), X(t) 


i 


Sunetion of 3 
relative velocity between equipment mass and 


aheii, in complex frequency aud time domain, 


respectively 
25 = geueral aotation for impedance function 
] = ghock wave decay constant, sec™? 
& = unotetion for characteristic equation 
Y" = gomua funetion 
Cy = Neumann Factor; 1 for n= 0; 2 forn>oO 
v = Poissoa's ratio 
X(T) = angle of shock wave envelopment of cylinder 
= cos” (1 - T) 
0 «/ mass density of mecivu 
A. = msgs density of shell material 
‘ i : : r Lpr “a 
‘© ach Clreular frequency of equipmen 
foundation model 
. = non-dimensional stiffmess of shell = a 72 / ac* 
O, = modal circular frequency of shell in-vacuo 
te = ghell stiffness in-vacuo 
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potential funetion of diffracted and radiated 


waves in fluid 
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L. INTRODUCTION 


A. Shock-Resistant Equipment Design 


in recent years routine tests have been conducted by the 
Nevy to measure the resistance of naval ships to underwater shock 
waves from non-contact underwater explosions. These tests have re- 
vealed that incapacitating damage to shipbcerd eculpments has occurred 
at shock intensities of only a fraction of those reguired to produce 
lethal huil demage ‘=! ) - At the present time, considerable attention 
is being devoted by various design anelysis methods to the predicticn 
of the responses that can be expected by equipments to typical motica 
iuputs. Because of the complicated interaction problems in the 
general case of a shock wave-shell-ecuipment interaction, shock motion 
inputs have been arrived at only after evaluation of vast amounts of 
experimental deta from reslistic underwater explosion tests. Retional 
utilization of experimental date requires correlation between test 
charge intensity, shock motion recorded, and shell-equipment inter- 
action parameters to include reletive mass, on-board location, and 
natural freqiencies. By realistic modeling of equipments for which 
dynamic response data from underwater explcsion tests are available, 
it is possible te solve for the shock motion input at some reference 
point, i.¢., inboerd side of sheli. A wide range of dynamic situations 
have to be investigated before a general set of input curves can be 
arrived at. 

An elternate approach is to tackle the shelil-shock wave 
interaction, making sufficient assumptions to make the problem 


solvable and to investigste: 
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(1) the qualitative validity of the solution 
(2) possibility of the imtroduction of experimental 
iufermation to adjust the qentit tive solution. 
Tais thesis shall investigate the reasibillty of senelytically 
describing the shock weve - hull - eguinpment interaction for the 
case of a submerine-type model. This description will take the 
form of en impedence. Puce to non-rigorous, yet physically meaningful, 


assumptions, the result will be termed a quasi~impedsnce. 


B. Shock Desisn Inputs to Several Analysis Methods. 

A design input may be defined as a description of the 
gnucek fereing or motion function to the design system. Dependent on 
the sygtem chosen, it could be a velocity or acceleration as a 
function of time at some reference point in the ship or some pressure- 
time description of the shockwave itself. A summary of the design 
inputs to the three best-known design methods are outlined below. 
1. Static Shock Design Number 

The most rudimentary of the methods, the input is a 
shoek design number, dependent upom equipment size, ship type and 


the attack direct ton") : 


It is applied to the equipment weisht to 
produce en equivalent force which the foundation and mounting should 
withstand. <A proposed modification to the approech would modify 

the number by a resiliency factor based upon the ratio of equipment 
natural fremency to the dominant frequency of the invut forcing 
funetion, as determined from peaks in the frecuency spectrum of the 


input motion“) ; 
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2. Wormal Mode. 

The normal mode approech currently deseribes an equivalent 
shock motion input (frequency domain) te a mathematical model at 
some inboard Tixed-base reference polnat, i.e., for a submarine, the 
inboard side of the pressure hull. The fundemental deseriptioa of the 
input is called the shock spectrum. The shock spectrum is defined 
as the trace of the maximum absolute response of a single degree 
of freedom system to the applied rigid base motion as the frequency 
of the oscillator is varied from 0 - ©. ts drawback as a design 
input for naval shock lies in the fact that it fails to account fer 
structural interaction between the coclllator and the foundation 
exciting it. The spectrum provides no information with regard to 
time histery, but rather an upper bound on the maximum value of the 
respouse ag & function of frequency. Extension of the non-carmped 


(7) and discrete-mass 


wg 63! 1038) 


vibration absorber from classical vibrations 
laboratozy experiments have led some obsexve to propose 
that the values of interest in a shock spectrum tend to lie in the 
Bulls of the plot, rather than at the peaks. A fiducial envelope of 
tae peeks is valid as an input when dynamic reaction of the driven 
component is mnable to affect the motion of the driving component. 
So as to account for the interection, the fundamental shock spectrum 
is modified by equipment weight, on-board location, and natural 
frequency. The result is termed a "design shock spectrum” and is 


the desiga input) : 


By assembling enough experimental data 
(equipment response) ranging over 2 sufficient band of equipment, 


natural frequencies, relative mass rabics between equipment-foundeticn 


wlie 


and ball mass, on-board location, charge (skock-wave) intensities, 
end ship types, degign inputs can be arrived at. Hach point on a 
Gesign input chart requires equipment-founcation modeling amd inverse 
normal mode solution for the foreing function required to produce 

the measured response. For the submarine case, a modification to 
this approach has been suggested which moves the reference point 
outside the sheil and adds hull flexibility to foundaticn flexibility. 


Hall resiliency would be experimentally erromaned ©? (43.) ‘ 


3. Iteration. 

Essentially a numerical method, it describes the bull - 
foundation - equipment interaction by a non-linear two-mass system. 
ft is solved by convergent iteration on the equations of motion. 

The design input is a pressure~-time history of the free-water shock 
wave which is the forcing function on 2 discrete mass describing the 


36 
hull “effective mage” 2°) ‘ 


C. Mechanical Imp edance Approach. 


Concurrent with renewed interest Im the sdaptaticn of 
classical dynamics to complicated structures, has been the adoption 
of en anelytic tool from electrical engincering known as the 
impedance techniqe. 

The use of impedance permits combinations of @ Large cor- 
plicated arrangement cf components to be hendled as e single entity. 
Computetions] efficiency is Increased because only input end cutput 

ata are handied, eliminating explicit introduction of ell the 
parameters characterizing each component "black box". Further, the 
parameters contained in each "bleck box” are 2a funetion only of thet 


particular component and are not influenced by the component's 



















aie 


(A PMEENEE (Ie ke) eee mre: ot neem Laat be 


Ke 5) OR ad a cet cnt cents emma 


Cate TU gl cate +51 ORR Sem ieee anton, & _—. all & 
et a 
} PARTS 5 ee sent ath ss caters ar) am ; 


COE CUNY Sl sere ice Mtneie eons ont uw ev 
wey. FO A WAM at i das de a 
, | 2 Linton i Abn cena oat 

Ht ah cme! Wiley Contemeey ested oe 
TE A eae a eae Ube a" 
he eh LT tite 
A NE At emai waterside) © 44 vey sid wa 
ieee Lo ee | 


= 





1 Ate ster AF nth Senet Eby iavicamay 
A et ii vette ies Soe lartuae Bl wy 


~— mee 





surroundings. The use of such e& technique in the analysis of shock 
motions appears ideal in that it provides: 
(1) a flexible procedure allowing for the introduction 
of empirical and quasi-analytic functions describing 
“plack box" interactions. 

{2) direct utilization of experimental data 

(3) maximum ease of evaluation of effects on the system 
response due to component changes, i.e., removal, 
alteration, replecement. 

To characterize the interaction of a “black box” with its 
surroundings, the impedance technique uses terms such as load or 
Griving point impedance, output impedance, and transfer impedance. 
A further refinement of the impedance techniqe formalizes these 
interaction impedances into two linear equations which relate the 
input point, 1, to the output point, 2, in terms of their four-pole 
paraneters: 

Fy = € Fy + @15 Vo 

Vy = @5) Fs + €55 Vo 
This approach is called the Four-Pole Parameter Technique and is an 


analytic tool in this investigation?) (23) 


D. Shock Wave - Hull Interaction 

Pressure wave interaction with cylindrical obstacles was 
duvestignted Classically wy Releigh’!*). the problem nas 
received revived interest in the application of the basic approach 
to obstacles whose thin cross sections are allowed to deform 


elastically. 
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The earliest of these studies ‘*7) used a modal analysis snd 
retained the full, exact form of the diffracted and radiated wave. 
The analysis was handied by transform methods. Due to the complexity 
of the solution (presence of complex poles and branch points), inversion 
by contour integration for more than several modes was tedious. Mindlin 
and Bleich’=”) py introducing an approximation valid for small tine, 
simplified the solution of MT) and presented numerical results for 
the first three modes of cylinder response. The simplification 
emounted to giving the propagation cheracteristics ¢ eae plane wave to 
the radiated and diffracted wave. Extensional effects in all modes, 
but the dilatational mode, were neglected. Baron ‘29) extended the 
analysis of (28) 


in all modes, using the plane wave approximation. More recently, 


to include both inextensional and extensional effects 


Haywooa@°) has contributed a derivation allowing for & more accurate 
approximation to the characteristic impedance of a cylindrical wave 
allowing for an “after-flow” effect. The approximation reduces to 
that of 25) for the case where the modal “after-flow" constant 

reduces to zero. This treatment extends the time validity range to 
times of the order of shock wave transit time, while reducing accuracy 
for small values of time. 

These investigations have been restricted to the two- 
dimensional problem, i.e., response to long shock waves. Longitudinal 
effects in the absence of cross sectional deformation are investigated 
in (33) ona") For this “whipping” moje, the model is a cylinder with 
rigid cross section, which is allowed to bend elastically about its 
longitudinal axis in response to a local (spherical) step wave. For 


the case of a symmetrical cylinder-wave situation (plane wave - cylinder), 
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the results of Murray compare in general with the translational 


(3 = 1) moge of Mindlin - aieacn 2°) , 


E. in-Vacuo Sheli Normal Modes. 

Tae natural frequencies for a circular cross-section 
in-vacuo, es applied in each of the previcus approaches, are 
derived utilizing LaGrange's equations by Timoshenko =) . TRe 
in-vacuo frequencies of vibration of a simply supported thin cylindrical 
shell, reinforced by equally-spaced, circular ring stiffeners, are 
derived by Gelietiy ‘3? : Bleich”) has derived an approximate 
method for calculating the normal modes of ring-stiffened cylinders, 
found to agree with tue more exact treatment or (3) within + 10 per 
cent. These methods satisfactorily describe the natural frequencies 
of a more realistic model for lower modes without undue labor. 


The general form of the results of (45) are listed in Appendix D. 
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Il. ‘THEORETICAL PROCEDURE 


4. Description of Systen Model. 

The hull is characterized by a simple cylindrical shell 
with a circuler cross section immersed in an infinite acoustic 
fluid. The response of the shell will. be elastic. Mounted at 
Point A on the sheli is a simple oscillator representative of a 
foundation - eqilpment chain, whose elastic response is to be 
determined. The equipment response is defined to be one-dimensional. 
Although the physical situation is defined to be two-dimensional 
(Fig. 1), longitudinal effects are compensated for in part by the 
uge of an empirical factor, called the effective length ratio. 

No "whippine” effects are considered. The ecuipment - foundation 
response time, while remaining a variabie in the solution qasi- 
impedance, is stated to be of the order of the shock wave transit 
time (5-6 msec.). The shell-simple oscillator combination is 
subjected to en acoustic, plane, exponentially-decaying step wave, 
normally incident to the longitudinal axis of the cylinder. The 
use of linearized theory is valid only for pressure waves which 
are not too intense and boundary excursions which are relatively 
amall. Syrmmetric modal deformation of the cylindrical shell. is 


assumed. 


B. General Approach. 
fhe sequence of events in the derivation of the cylinder 
quasi-inpedance will be outiineé before starting the muthematical 


details. 
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Define the shell response, the system input force 
and velocity, and a potentiel function for the 
radiated wave in terms of Fourler expansions of their 
generalized cocrdinates, generalised force and 
velocity, and modal potentlel fuacticn. 

Define en effective modal. mass of the shell. 

Setup of the equation of motion of e& point on the 
shel. The shell. is rsea@ in the median and the 
pont feels the Interaction with sn equipment - 
foundation chain mountec inside. The omeli is acted 
upon by the summation of the inacident, radiated 

ani diffracted pressure waves. 

Description of the naturel moGel frequencies end 
shapes in-vacuo for the shell model chosen. 
Approximate the characteristic Impedance of the 
radiated wave. 

Prescribe necessary Doundaryy conditions. 

Using (5) and (6), eliminate the fluid potential 
function from the medal equation of moticn. 

Express the modal equation of motion in terms of its 
Laplace trensformms and convert to normal four-pole 
pareneter Toxmat. 

Set up the fourepole parameters for an appropriate 
model of the foundation - equipment chaln. 

Using mitrix techniqies, combine the shell end the 


equipment - foundation impedance functions. 
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ll. Invert the resulting Laplecian expression of (10) 
to the time domain sclving for the equipment re- 
sponse to one mode of the shell response to the input 
shock wave. 

12. Sum the responses of the equipment to the modal 
excitations of the shell for as many modes as 
required. 


13. Convert to response relative to hull if desired. 


C. Impedance Parameters. 
A general discussion of the mechanics of forming cuadripole 


parameters for a wide variety of paysical systems is available in'®3) ‘ 


The brief discussion included here is soleiy for the sake of continuity. 


1. Definitions 

A linear system which possesses a single input point, i, 
and ea single cutput point, 2, may be described by a pair of linear 
equations of the form 


F. + €,, V 


FL * &y %a* “ye "2 (1) 


V, = &5) Fo * €n9 Vo 
where €117 "12? Sop and €,, are the quadripole or four-pole 
parameters of the system. Solving Equation (1) for the four poles 


gives: 


0, i.e., stetion 2 blocked; 


% 

i 
sh 
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= 0, i.e., station 2 free; 
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V 
€ -= with V. = 0, i.e., station 2 blocked; 
2 


lL 2 
"7 
G55 = ¥, with Fy = O, i.e., station 2 free. 
The 24 parameter is thus the force trensmissibility from station 1 
to station 2 where station 2 is blocked. Parameter E15 is the 


transfer impedance between station 1 and station 2 with station 2 
free. Parameter fn) is determined by the transfer impedance 
between stations 1 and 2 with station 2 blocked. Parameter Eno is 
defined in terns of the velocity transmittance from station 1 to 
station 2 with station 2 free. Eoustions (1) are conveniently 
expressed in matrix notation as: 

a e — ®12 ad (2) 
7 Se Coa ‘e . 
The first bracketed term on the right is called the cuadripole 
parameter matrix. Matrix algebra provides a formalized and 
precticel method of forming compenent quadripoles from fundamental 
elements, and of obtaining the overall system quadripole from its 


components. 


2. Quadripole Impedances. 

The elements of the four-pole matrix are naturally 
related to the various types of mechanical impedance that describe 
the component. The relationship between the two are listed below 
in order to provide correlation between the two methods of 
nomenclature. 

The load or driving point impedance characterizes the 
elements that follow iit: 
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is dependent on both the quadripole and on its load impedance. It 


is related to the load and the four-pole by 





The input transfer impedance is 
ts i = 2 Z, & 
219 ° "a S* “a” 


The force transfer function is 


EF 
wh = + €,,/Z, - 


The velocity transfer funeticn is 

v,/V, = Gan * Gp, 2 
3. Forming the Parameters. 

For some simple models direct combination of the know 
parameters for fundamental elements (Appendix C) using matrix 
aigebra will produce the quadripole for the system. The definitions 
given above also provide a means of introducing experimental data 
into the analytic solution. However, in the general case it is 
necessary to: 

(1) define the system 

{2) describe the system in teras of a single pair 


of input end output junctions 
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(3) sclve the performance equation using operational 
calculus methods 
(4) vesolve the equations into the form of equation (1) 


ex (2). 


ent Interaction 


D. Ausalysis of ShelleShoeck Wave - Foundation - Equi 
a. Incident Wave. 





The input is defined as an exponentially decaying step 


wave of the form 


P, = P. ent u(t). 


1 
The simple step is translated across the cylinder with the use of 


a time variant delay factor 
si = ( 1 - cos @) (origin taken at cylinder axis) 


g0 that 


? enlt = 7) 


= P, U (t - T) (3) 


i 
The total pressure and the radial component of the particle velocity 


are expanded into Fourier series of the forn 


& 


P, =) Pao cos n 8 
n=O 


€d 


U,=) U,, cos ne (24) 
n=0 


The Fourier modal coefficient P,. is obtained from equation (4) by 


in 
modal orthogonality after miltiplying right and left sides of 
equation (4) by cos n 8, integrating from 0 to 2x, interchanging 
the order of summation and integration. The integrand is zero 


from © to 2n by the definition of the step, u(t - +). 
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This procedure yields 


3(%) 
P, (1,%) cos am & O< & (t) < a 


3 tj 
© S-~-5 


where @ (t) is the angular position of the incident wave. Using 
the Laplace time - lag theorem for translation in the time donsin, 
substituting from equation (4) and transforming with respect to the 
non-dimensionel time Tf = et/a yields, 








e = eos (1-T) 
Poy b fery [oT eno 029 <a 
(5) 
Similar development appiies to the radial modal particle velocity. 
Introducing an additional cos 6 because of its radial definition 
and using the characteristic Impedance of a plane weve, 
6 (tb) 
1 41 ' PA 
Uy, * +3 J le (2,%) cos n 0 eos 943 O€ 9 SB 
0 
and 
9 (t) 
(s+ By, Saxe ~ fle" ® 5 ° cog n 0) cos 9 a9 OS 5S & 
0 (6) 


At @ = 2n, the integral in equation (5) becomes a definition of 

the modified Bessel function. During treusit time, they are 
referred to as incomplete modified Bessel functions. Their numerical 
evaluation is discussed in Appendix B. Using this definition, 
equations (5) and (6) reduce to: 
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(prime denotes derivative with 
respect to zg.) 


2. Interaction Pressure Field. 

Contact between the incident wave and the cylinder pro- 
daces scattered and radiated acoustic waves whose potential 
function must setisfy the wave equation. The scattered wave 
arises because of the presence of the cylinder. The radiated wave 
iS cansed by the response motion of the cylinder. The scattered 
and radiated waves are assumed to move radially ont from the 
cylinder. The radiated pressure and velocity are defined in terns 
of the potential function % by 


5 
Pr =P 3 


as 
er” ” Or 


Modal contributions, corresponding to the in-vacuo 
generalized coordinates of cylinder normal modes, are obtained 
by expending the velocity potential, this P. and Uy. into & Fourier 
cosine series in 6 of the form of equation (4). Equating coefficients 


ef the expansicns, yields: 
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wnere the generalized potential functioa also satisfies the wave 
equation in cylindrical coordinates 
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3. Cylindrical Wave Approximation. 

The use of the complete form of the weve equation above 
to provide a relation between P_ and U_ (characteristic impedance) 
yields a solution whose inversion fram the complex freqency domain 
is prohibitively tedious. 

An epproximation by nagwoua' 20) to the characteristic 
impedence of a cylindrical diverging wave yields results valid 
for times of the order of the shock wave transit tine. Heywood's 


derivation yields: 


ar cas hCUF Ch 


or 


The modal constant g, is the average function whose limits, while 
time ranges between 0 and ~, is: 


O<@ < 1/2 


1/2<¢e,<2 forn >i (9) 


The main effect of the approximation is to give a fimal veicelty, 

or “afterfiow" to the fluid, following the passage cf the pressure 
wave, whieh goes to zero only after the vreegure has fallen below 

its equilibrium value. For distances infinitely far from the 


cylinder, equation (8) reduces to the plame wove approximation. 
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kh, MNetural Frecuencies of Shell Model. 

For survey purpeses, the simple circular unstiffened 
shell mo@el is considered adecuste. Extensional effects in all 
modes except the dilatational (a = 0) are neglected. The 
derivation?) (15 dor the nozmmei mode frecuencies for a simple 


eireulsr cross section yields: 








a," rd E = (1a = 0) (10) 
Pe i 
with § = ——, 
(ley) 


2 2 = 
2 1 i T 
i (inextensional - flexural) = ati =a) x 5. x z x mp 
lin a o 


The natural frequency of the translational mode n = 1 is, of course, 
equal to zero. 
In order to reduce the order of the characteristic 


equation for a later approximation, we now define 


us 
= os as the model shell stifmess, 


where Ww, = m, (n= 0) and 


@ 
= Gp) a (m > 0). 
a 





ad 


5. Equation of Moticn at Point on Shell. 
Summation of the modal forces acting at point A in 


Fig. i gives: 


°% .5*—-8 + Bat Py (12) 
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where q is a generalized coordinate and is defined in terms of 


the total radial deflection w by: 


@ 
we) q, cos a 6 


n=0 


6. Boundary Condition. 

At any point on the surface of the shell, the modal velocity 
of the shell mist be equal to the net particle velocity of the in- 
cident and radiated - diffracted modal waves 


4 * Van > ow (13) 


~. Net Shell Response Function. 

An operational flow chart combining the preceding steps 
is given in Figure 2. 

Transforming equations (8), (12) and (13) by Laplace 
trensforms with respect to time ylelds 


= 
ms (t) 


a ) B= Fay - Boy + See (yy - &) iad 





Rearranging and changing the time variable to a non-dimensional 
time T (= ct/a) yields: 


) A= Fe [oles BY Py oo Bay 
+ s* pe 0. (15) 


where 4, the characteristic equation of the shell response, equals: 


e+ GRE rans Sh) o + SG, (16) 
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From equation (7), 


therefore, equation (15) reduces to: 


pe (A) a =z «8 (s + é.,) Fy + | a(6 + g,) + 22) Pan 
"Qn 
Introducing the definition of Pen for step decaying shock wave 
equation (7), we obtain 


. s (6 + B)(s + g,) 
. (els (s + @) 4, (s)+ s° 1. (s) | = 


“ ee a = Q (18) 
's (s + é,) Too {(s) + a” .. (s) je n 

E. Equipment - Foundation and Shell Four-Pole Parameters. 

In general, there is no limit to the complexity of the 
model of the equipment - foundation chain mounted at A in Figure l, 
except that it be capable of representation by point attachment. 
For the purpose of this investigation, a simple oscillator was 
chosen. The mass and natural frequency may be varied to aeveruine 
thelr effects on the shell - equipment interaction. 
The mass four-pole ag derived in Apvendix C is: 

Bo] 


The spring four-pole is found to be: 


lS/k 4 


They are connected ia cascade, thus are combined by matrix 


mltiplication to give: 
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This four-pole has been transformed with respect to time. In 
order to conform to equation (18), we change the time variable 
to non-dimensional form T (= <), The four-pole decones: 


Cc 
2 2,2 
P Ms". e/a 
2 1+ 


(20) 


% $e 


We see that equation (18) fits the general form of one 
of the linear four-pole equations (1). There are several points 
to be discussed or rationalized before we give (18) an impedance 
tag. 

(1) Equation (15) deseribes the modal response of a 
distributed parameter system to the spatially 
aistributed, time variant force of the pressure 
loading and a dimensionally siniler equipment reaction 
force. Implicit in the impedance concept is the 
necessity for point input - output description of 
the system. The application of the impedance 
concept can be rationalized if we consider as the 
imaput to point A en equivalent poin® force which 
would produce the same response at A as the distributed 
pressure loading. Attached to the output of point A 
is the equipment model. 

(2) To further satisfy our point requirement for the 


eross section, we stave that: 
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(4) 


(5) 


P 
fh, (equation (15) ) = = (21) 
where ‘i. is an eqivaleat point force for the 
cross section, but remains a distributed force 
longitudinally representing the equipment 
reaction foree per undt of longitudinal length. 


K = arab (22) 


efT 


where L. £ is defined to be the ratio of the 


f 
longitudinal length cf shell feeling the influence 
of the equipment reaction, to the actual equipment 
iength. Love 48, Of course, an empirical input. 

It will be a function of relative stiffnesses and 
masses, but a range of 2-5 1s considered reasonable. 
The interaction function is defined here only for 
an exponentially decaying step wave. However, 
physically, this is qite acceptable, since by 
superposition and the use of the Laplace delay 
factor, the characteristics of any nermal pressure 
aouGing tan be silmilaved. 

The net impedance function is defined in terms of 
en intinite series of modal impedances. 

The complexity of the model prevents complete 
description of the four-pole parometers. The 
problem is in essence siviler to the transmission 
line.problem in electriesi engineering. There, 
plane waves travelling along a trenesmission line, 


described by a characteristic ilupedamce, meet an 
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iaterface or loed impedance different from that of 
the transmission line. Due to the interaction, 
stending waves are set up in the trenszission line 
altering the voltage - current relatioaship (input 
impedance). The analogy is made by regarding the 
transmission lime as the medium surrounding the 
cylinder, and replacing voltage and current, by 
ferce and velocity. The essential difference and 
reason that the analogy cau be extended no further 
arises from the fact thet with the cylinder the 
problem is no longer bileteral. The input wave 


Impedance differs from that of the cutput. 


Iu order to form the four-pole parameters according to 


the standard recipe, we would first define the system components 


end component point inputs and outputs as: 


b) 


outboard 


FIGURE 3 
SYSTEM FLOW CHART 


Plane wave characteristic impedance = ec 


~—= = $C COS *£ 
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Input transfer imedence for the water component is: 
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ce. Foree transfer impedence for the water is: 


ad. Velocliy transfer impedance fur the water is: 
Uy 


eerie 
U; u,. 


e. Load impedeace at the water output junction is: 
Fo, 


a 


ce 


Lt is possibile to describe 





P P 4 
Pe? ad oy ~ UY 
or 
i. ue an 
Ps + Ye U 4 vy 


using equations (18), (7), (13), but we have not adequately described 
the interaction to be eble to define both without obtaining one 
redundant yet. However, for cases where the response of the last 
element in a chain is desired, i.e., output end is free, force equal 
to zero, one of the two four-pole equations is edequate to describe 
the motion response. All interactions enroute are stiil accounted 
for. Such is the cage in this examplc. 

Thus, we have not described the system components entirely. 
In fact, we have shunted the complicated lateraction at the shell 
surfeee and have dealt with a system whose limits are defined at 


the Input by the free stream shock wave and at the output by the 
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inboard junetion to the equipment - foundation four-pole. 

With these rather restrictive conditicns governing our 
transfer function, we call it a quasi-impedance. 

The shell function equation (18) is now combined with 
the fourspole of (20) and definition equation (21) to preduce after 


some rearrangement: 


Se. Py [(s + &) I. (s) +9 It. (s) le“® 
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a a) EG RE Ewe aA 
(23) 
xn 
- Bhi cs 2, 8.52 
where Beet Geet h)s + (“3 G,") 6 + 5H, 6. 
pc e 
(24) 
Finally, after defining 
+. * = = mataral frequency of (25) 
a foundation - equipment 
°.. FP 
V = % Apes a exis _— , 26 ) 
a" e 1 
“= *_ * 8, “5 x =, 
a? a” © ach ' 
nd expending 


7. (s) = ° v, | (s + &,) 1, I (8) tal. (s) le" 


i Nae a en 
a + B)(s° +o +8 3° + 28° +asta ) 


3 
(27) 
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dy Oy ii 
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and the reots of the fifth order characteristic equation equal to 
S59 i = 1 3 5 
After eliminating the Bessel function derivative in 


(27) using the following recurrence relation: 


It (as) = $5 I, (os) =a 0, (08) - 24, (cs) — (28) 


the velocity is inverted from the complex frequency domain. Use 
of the residue theorem from complex variable theory (4) makes the 
procedure quite straightforward. 
The velocity response for all modes other than translational 
(n = 1) takes the form: 
Pe 1-T 
a (8) [Gq = 2 = 8) yp (8) - (8) Ty 1){-2) BO) 
V. (2) = 
3 ao B){s, + B)(s, + B)(s), “ B)(s, + B) 


1. 
tn Yo 


4=1 


cars - T) - 


where H, = (s, v B)(s, - s,)(s, - 3, )(3) ” 8, (3, - 8.) 
H, = (a, + B)(s, - 8, )(s, - $,){8, - 8,48 = 8.) 


H, = (s, + B)(s, - 9, )(s, - s,)(s, - 8, )(s, - 85) 
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Hy, = (3), 4 3)(s), > 3, )(s), - 5, )(si, ~ 8 )(3, - S-) 

i, = (6, + B)(s, ~ s, (s, - 85 )(3,, ~ $,)(s, - 8) ) 

In the treaslational mode (nm = 1), the shell stiffness, 
a goes to zero. Tae coefficients of the characteristic equation, 


2743 2 
8 (s + &), 8 *@, 8+ @,), 





reduce to 
a he 
5°55 hat = (30) 
- sg i “ach 8 ¢ 
®2 = “3 (nach Sa * 98 “Ky ia, ) 
a, = 0 
Gy) = 0 


The calculation of the residues mist be repeated die to the double 
pole at the origin. The translational velocity of the oscillator 


mass is foun@ to be of the fora: 


(s, = 2) Ty. (0) 
Van > rev, [yt ]a- ¢) 
2 
62y, “ya ae 


1 + BYs, + Bs, + B) 


Foy, (let a= 28s Tha (Oy) * 81° Foe (6) os, - 0 
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where H, = (a, - 8,)(s, - 8.) 
A, = (5, . 8, (3, - 8.) 


Hl, = (8, - 63)(a, = 83) 
The velocity of the equipment mass es derived here is in an absolute 
reference space. The absolute velocity is 2 satisfactory measure 
of the effects of system and component parameter variation. However, 
from the stendpoint cf equipment damage criterion, a more useful 
quantity is the relative motion between the Ill and the equimment. 
Relative motion will, of course, gevern the maximum stress, deflection, 
and acceleration that the equipment ang foundation will see. The 
power of the impedance technique is revealed in the simplicity with 
which the relative velocity is calculated. 


(1) Equation (27) gives the equiyment modal velocity 


in frequency (s) domain, 


Pp 
° 
Von (3) = rac) 
(2) Equipment - foundation four-pole parameters given as: 


‘. P %, ye Pan - 
‘ 2 
Von ‘ “we 3n 

(3) Since Foy = Q, 
Vo, (8) = 55 (s) V5, (s) (32) 


(4) Therefore, 


P 
Tt; e ;' » ae. 
Von (s) Cao (s) x Sa) 
(5) Inversion to time domain, if desired, is identical 


in procedure to that of Van, (s). The poles of the 
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response function are of course identical since 
the characteristic equation for the system is 
unekanged.. 

(5) The relative velocity is then: 


& {s) = Von (s) - Von (s) or, 
x (2) = Vn (Tt) - Von, {f). (33) 


F. Numerical Procedure. 

Equations (29) end (31) will form the basis of the time 
respoase calculations. while the entire procedure is straightforward, 
the number of cycles required to trece the time history, coupled 
vith the superposition of an indeterminant number of shell modes, 
and the large oemounts of complex number arithmetic, lead to the use 
of a digital computer. Since the general form of the residues is 
determined, the machine needs only to: 

(1) aceurately determine the roots, real aad complex, 

of the characteristic equation. 

(2) perform numericsl quadrature or series approximation 
(Appendix B) to the incomplete Bessel functions at 
steps sufficiently small to provide time history 
desired 

(3) evaluate equations (29) and (31) as functions of 

time 


(4) superpose modal responses. 


fhe frequency response (s = iw) of the quasi-impedance 


is readily available by choosing the range and discrimination necessary . 
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Cuce satisfectorily programmec, ail parameters can, of 
course, be varied, using systematic procedures to determine 
existing trends. For the initial computer runs, typical parameters 
were arbitrarily chosen to characterize the shell model, and skock 
weve. Three senarate combinations of equipment - foundation model 
rass and noturel frequency were chosen in a typical range to order 
to evaluate the ability of the quasi-impedence to deseribe intersction 
effects. 

So az to be able to check by hand initial calculetions 
using & closed form rootefinder, the order of the characteristic 
equation was reduced from quintic to cquertic. This reduction was 
accomplished by neglecting the shell mass. The maxiwm error will 
be experienced for times shortly after the ghock wave arrival. 

The natural periods of the equipment models chosen, are of the 
order of transit tine. However, it is obvious from equation (23) 
that the camplete interaction is a function of relative hill - 
equipment masses. Therefore, no physical justification is offered, 
but rether the argument that the general ability or imeability cf 
the function to characterize the interaction will be retained. 
Details of the programning procedures are listed in Appendix A. 


The computer lenguage is IBM FORTRAN and the computer IBM 7090. 
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IIL. DISCUSSION OF RESULTS 


TE 


A quasi-impedance function has been derived for e« simple, 
but physiceily reascnable model. The function vossesses flexibility 
for extension to analytic models of a more complex nature without 
changes in its basic forn. I% further provides an opportunity 
for incorporation of empirical data in the form cf measured 
component impedances end shell geometry varieabies into an analytic 
description on a rationzl basis. The function describes the 
equipment model response in terms of nhysically meaningful paraneters: 

@. equipment - foundation moss, stiffness, and natural 

frequemsy or experimentally determined impedances 

». shell mass 

ec. shell in-vacuo natural frequencics 

d. shell geonetry 

1. radius 

2. cross-section moment of inertia 

3. shell thickness 

4, frame spacing, area end moment of inertia 

with the use cf Appendix D 

5. empirical input of “effective length ratio” 
e. shock weve characteristics 

1. peak pressure 

2. exponential decay constans. 
Complete evaluation can, of course, be made only after 


completion of 9 program of systemeatie parancter variation in 
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conjunction with comparison with measured field data from realistic 
tests. The munerical procedure to test the general validity of 

this approach during this exploratory phase is outlined in II-F 

end in Appendix A. Couclusion of the time allotted for the numerical 
example occurred before the digital computer program was completely 
Functional. Complete reliability of the root-finder subprogram to 
obtain the real and complex roots of the characteristic equation 

was not achieved. Two methods using the digital program and another 
hand-calculated test method were used in searching for the roots. 
Comparison of the results from these three methods lead the author 
to conclude that s combination of absolute end relative root 
magnitudes in the lower modes and relative root mgnitudes in the 
higher modes both prevent convergence and introduce mumerical 
inaccuracies into the root-finding technime. While the possibiilty 
of programmer procedural error is not precluded, careful checking 
has revealed no blunders (as opposed to errors). 

The majer program using a double precision arithmetic 
method for automatic commtation in a complex mumber mode of the 
pole residues and time variance appeared to function properly. 

The failure to achieve general program reliebility naturally 
prevents inclusion of the velocity responses calculated. However, 
for the modes whose chun tenistae equation roots were reasonably 
accurate, responses were of a damped oselllatory form as would be 
physically expected. As the weight and natural frequency of the 
equipment were varied, response calculated conformed as expected 


(25); 1.e., the larger the mass, the lower the peak mass velocity 
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and the larger the rise time. The influence of mass of equipment 
on peak velocity can be seen directly in equations (27) and (29). 
The direction of the mass effect on the rise time is less obvious, 
but it is clesr that the mass influence is directly proporticnel to 
its magnitude. The accuracy of ereluetion of the incomplete Bessel 
function for large negative root values during early times of the 
transit period became progressively worse as the mece order increased. 
The cause of this is shared by the massless shell used to simplifv 
roct-finding problems and the form of the equation for small 
times involving smell cifferences between large numbers. Accurate 
results for these early phases of the response motion can be 
achieved by series expaasion of the impedence in s end subsecuent 
tern by term integretion. The quasi-impedance would then be 
Geseribed hy a two-phase function, each phase with a time validity 
renge. The need for this added complexity needs to be evelueted. 

Equation (27) describes the velocity response relative to 
a set of fixed coordinates, which 1s adequate for scme parametric 
studies. However, in most situations the damage evaluation or 
*“verformence level” is better deseribed in terms of response 
motions relative to the hill. The utility cf the impedance approach 
is demonstrated by the simplicity with which the relative velocity 
is obtained, as shown by equations (32) end (33). 

In general, the investigator is corvincea that the system 
mode]. is ultimately capable of providing 2 qalitative and quantitative 
estimate of the response of a real cross section cver & range of 


test conditions. It provides en analytic model bearing some of the 
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characteristics of a distributed masa system with the convenience 
of a discrete impedance chain system. As such 1t provides a means 
of attaching quantitative meaning, in the form of time response, to: 

(1) skoek-spectrum avull seen at equipment Pixed-base 
natural frequency In discrete mags experiments 
lovestisating distributed mass phenomena. 

(2) knowledge that shock spectrum ~ normal mode 
approach, by destroying time history of input motion, 
results in over-conservative estimeted response (due 
to neglecting of moGal time to peak in surmation 
precess ). 

By describing the response in terms of physically meaningful 
parameters, the author believes thetthis approach will allow a 
better understanding of the basic phenomena. In providing a method 
of breeking down the complicated system interaction into component 
interactions capable of description by engineering approximatica, 


the impedance approach is cf great value. 
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CONSLUS IONS 


lL. It is possible to formulate an analytic approach describing 
the water applied pressure loading ~ fmlil - equipment inter- 
agtion with the use of engineering assumptions, capable of 
solution with a reasonable amount of effort. 

2. The use of mechanical impedance and transform methods provides 
a@ powerful analytic tool and a flexible format and procedure 
for the introduction of more exact model components and of 


expirical data. 








Systematic parameter variation of the derived quasi-impedance 
function should be conducted to provide rational basis for 
the establishment cf design inputs and goals. 
The impedence function should be expanded to include the 
effects of: 
&. multi-degree of freedom equipment - foundation model 
b. fovndation and equipment damping 
ec. longitudinal whipping response 
d. experimental - empirical inouts 
Response variation should be investigated as the equipment 
attachment point varies around the periphery. 
Investigate need for two-phase impedance function to extend 
accuracy during early times of transit. 
In conjunction with 2e and 3, investigate the effect of attack 
angle of incidence. 
A quantitative investigation of the “spectrum dip” phenomenon 
for a distributed mass model of the form developed Bere 


should be conducted. 
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APPENDIX A. 


DETAILS OF PROCEDURE 


1. Computer Procedure. 
8. In crder to expedite hend calculation checking of the 
IBM 7090 computer pregram Guring the exploratory 
phase of the problem, the order of the characteristic 
equation A is reduced by neglecting the shell mss. 


As a resuit, A reduces to: 
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The roots to the charecteristic equation above will 


be called fy Sp 349 3), 


Invergicn of (A-1) ia accomplished by calculating 
the resicues about each of the poles, s i: 
The generel form of the time response for ail modes 


other than translational, is: 
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(2 -a+s, )Z,. ( * 
p {) nee a = Roanmieney ‘le ~83 (1-7) 


(a,-m+s,j1,. (85) + “2 a Oe (s,) jee) 
bd 81) 


* V5 Gy (89) Get BG IE 
+, ty be, ete ats. a Stee el _-83(1-7) 


+ 


(g,-0*s; i, (3),) + 5), 1 = (3; ) - 1- 
Vo Sq (8) ent + B)(s, = 3, )(s, = =) * za oni 


(ma x 1.) (A-2) 





db. Shell and eqiipment - forndation parameters were 
introduced as data to the machine, which calculated 
the modal coefficients of the characteristic equation 
of {A-1) for the first twenty modes. Subroutine 
rootefinder was progranmmec: to seareh for the real - 
complex roots of the resulting quartic. The first sub- 

rogrem used was en IBM SHARE Distribution Progran 

#1124, “MULLER"*. Persgisting machine overflows, due 
perhaps to coefficient incompatibility, resulted ia 


the author prograaming the Euicr - Descartes Method 


*"MULLER” is based upon the theory of Miller, J., “A Methed for the 
Solution of Equations with Autumatic Computers”, Mathematical Tables 
of fztomatic Computation, 1956, pp 208, 225. 
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for ciosed form solution of quartic, previously 
used to check computer results. The basic steps are: 
1. reduce quartic 


ae ree Ee ee 


to form 


- + ax® + bx +e=0 
by. substitution 

y = (x ~ p/h). 
We find 


a = 29.5 5P 
3 
— l6ér - 8pq + 2p 


2 


- . s + 16qp"_- Ghrp + 256 s 
25 


2. With £2, m, and n as the roots of the 


resolvent cubic: 
2 2 
3 8, ,2 a = he mm < 
eee eee sae 


The reduced quartic roots are: 


x si(-/e -Jm -Ja ) 
X= + @/ip Ge + ) 
* "<5 


44 ‘a -/@ te ) 
mar ve - /m -J/a } 


where upper sign is used if b is positive 


and the lower sign if b is negative. 
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c. Having calculated the roots of the characteristic 
equation, the computer was programmed to evaluate 
the modal incomplete Bessel funetion qs of equation 
(B-1) at three-degree increments in envelopment angle 
with a three-station Simpson's Rule. After storing 
these values, the residues of equation (A-2) were 
evaluated at intervals of the non-dimensional time T 
corresponding to the envelopment angle (6) increments 


of the integrated Bessel function. 


2. Typical Shell end Equipment - Foundation Parameters Used. 


a. Using typical values for E, I, a; 9; P52 Wh» 


yielded for the simple shell, 








2 
w,? (extensional) = (2.13 x 10°) 28%, 
O 2 
sec 
2 2\e 2 
w (flexural ) -aGsay x 782 oe ; 
L+a SEC 


a? = 4.23 x 10° 1p,/ft* gurface/ft Length 


an = 3035. (1 - n*)* 1p, /#t* surface/ft length. 


bd. Three combinations were chosen for the equipment - 
foundation parameters: 
1, M, = 155.2 slugs/ft. 
k = §.515 x 10° 1b,/ft. 
f = 30 cycles per second (cps) 
2. M, = 7-76 Slugs/ft. 
k = 1.960 x 10° Lo -/ ft. 


f = 80 cps 
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3- M, = 155.2 sliuge/ft. 
k = 3.2x 10° Lio, / 2. 
fF = 120 eps 


Sheck Wave Characteristics. 


A peak pressure of 2000 psi and a decay constant 


1 


8 = 20.8 sec ~ were introduced as data to the 


machine. 
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APPENDIX B. 


NUMERICAL EVALJATION OF INCOMPLETE 
AND COMPLETE MODIFIED RBESSEL FUNCTION 
OF THE FIRST KIND 


During the shock wave transit across the cylinder, 
evaluation of the Fourier modsi coefficients of the incident pressure 
yields the function: 

@ 
Ino (s) =< | e® $8 © oo no ac (B-1) 
. Of ea 


where 6 ig the angle of envelopment. For times greater than the 


transit time of the shock wave, 15 becomes: 
ert 
L 6S cos 6 
‘ (s) = ¢ e cos nd do. (B-2) 
0 


Equation (C-2) is the definition of the modified Bessel function 
of order n. Equation (C-1) is therefore termed the incomplete 
modified Bessel function of order n. Equation (C-1) is not integrable 
in closed form, but is capable of evaluation by: 

(1) mumerical quadrature 

(2) rapidly convergent series approximation. 
Figure B-1 plots I, , (s) as function of Tf (or equivalently ©) for 
a specified real root of the characteristic equation, as obtained 
using three-station Simpson's Rule. It is presented to give an 
idea of the gress form of the function. Im general, 5s is complex 
and requires integration of real and imaginary parts. Equation (B-1) 


takes the form: 
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: (s) = | e (erly) =e s cas no a) , which can be 
0 


written by Euler's identity as: 


8 
= | sows | cos (y cos 8) + 1 sin (y cos 5) eos nd ao , 
e) 


or 
S) 
2 [ e* £98 © cos(y cos 6) cos no ad 
0 
5 
+4 | e* COS © oan (y cos @) cos no ae (B~3) 
Q 


Where x and y refer to the real end imaginary parts of s respectively. 
A digitel machine mekes the evaluation of (B-3) .1ite effortless. 

An alternate method cf evaluation is to use the Jacobisincer 
formula and a serles expansion of the modified Bessel function. 
Bessel functions of non-negative order n = 0, 1, 2.... are single 
valued integral functions of s. They may be formed as cocificients 
of the Fourier series 

=e 


7, = x y! cn i, (s) cos 2% (Jacobd-Ancer Formila) 
=0 (B~4) 


where ¢,, (Neumann factor) = 1 for n = 0 
= 2@forn=1, 2+». » 
Using (B~), (B-1) goes to: 
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The c., 1, (3) 43 not a function of the integration, so that we can 
bring it outside the summation sign and interchange the integration 


and summation signs yielding a form: 


) ¢ 
TT Ti9'8) T3(s) | cos nO do +2 I, (s) | cos no cos 6 49 
6) 0 
§ 
+ al, (3) | cos n9 cos 26 a6 
0 
6 
< al, (s) | cos no cos 35 ao 
0 
+ esase rs (B-5) 


Equation (B-5) involves for real s only tabulated functions and 
integrations soluable in closed form. Further, because of rapid 
relative decay of the modified function with increasing order, 
five or six terms should be adecuate for engincering accuracy. 

Yor non-tabulated arguments the modified Bessel 
function of the first kind, 


nt 
I (s) =i J (is), 
where J (s) ig the Bessel function of the first kind is defined 


by: 


1, (s) -) (x/2)* ** (B-6) 
a (xk + m)! 


gm 2 ry 
gt r (m+1) 2(2m+2 2.4#(om-+2 jl amt 
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with © (n) gamma function = (n-1)! = (n-1) © (a-1). 
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For gmall values of s, (B-S) cam be approximated by the leading 


terms: 


— 
i, (9) + o (B-7) 


For large values of s, 


5 
(2n s 


In manipulating the modified Bessel functions, the following re- 





currence relations are useful, 


= I, {as)=al, (as) - ~ I, (as) (B-9) 

Sot (as)=at., (as) +21 (as) (B-10) 
as “m m1 5 ‘ 

ei fe s)=al[ (as)+2Z.. (as). (B~22) 
ds “m * a m+L a 


rl (a s) ~ rel (a s) = £ I. (as). (B-12) 
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APPENDIX C. 


FOUR~POLE PARAMETERS 


lL. Development of Four-Pole Parameters for Fundamental Elements. 


&. Mass 





Performance Equations are: 


Ve ¥,-¥ 


FL - Fi, = sal 


where F and V are arbitrary complex functions of 
time and F and ¥ are the transforms of these functions. 


It is clear that the four-pole matrix is then: 









‘lL om 
fm} - 3 7 
b. Massless Spring Four Poles 
Fy } Vy Fa, Yo 
2 


Perforasnuce Equetions are: 


Fe Fi, = Fo 


k 
PF, =k | a xp | 3 (Vy - Yo), 


§ 
= oh + 
or V, = § Fo * Yo 


The four-nole matrix is then 


_ 4 
ik} = ia/k Ll. 


ec. Viscous Damper 
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Poxrtormencc Equations are: 


a. 
a + V 


i 2 2° 


The four-pole matrix is tes 


e) 
oD ae 4 


2. Combinaticns of Camponent Quadripoles 


a» 





Series Combination. 
If the output of one component is connected to the 
input station cf another system, the two systems are 


in series. For example, if the two components ocicw 


are 


» 2 Sos F,, V 
®@ | 12! © 


comnected at point 2, the components are connected 


3 


in series. If the four-pole of fi} is 
A >| 
cD and 

the four-pole of {2} is 


“Be F 


q ql? wea the 


fourepole relating Fy and vy to Fs eud V3 is 
(A BIS FI 
le 4 r H! . | 
So that the net series connected system four-pole 
is equal to the matrix product of the four-pole 
matrices agsoclated with each of the systems making 


up the composite systen. 
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Parallel. Combination. 
Components are considered to be connected in peraliel 
when: 
(1) all the input and output junctions move 
with the same velocity 
(2) the net input/output force of the resolvent 
quadripole is equal to the sum of the input/ 
output ferees of the constituent four poles. 
If the system obeys Maxwell's Reciprocity Principle, 


then 








The matrix s of the composite system relates input 


to — by 


i] + o ~ 


wnere [5] = [I “Ais 
iD, /By “1/3, Pele ~1/B, 





[2/B, A, fe, * 1/B, A, /B, 
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Appiicaticn of the above formilas to a spring-damper 


parallel combimation yield an equivalent four-pole 


- Of 
S ‘ 
rok 4. 


More extensive developments are presented in (23), 


(1), ame (2k). 
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APPRMETX D. 
CPOE Be A ee 8 ET ETI 


DI-VACUS NORMAL MODUS OF RING-~SPIFFENED SHELL 


Analytic methods sre available providing az additional 
measure of realism in the description of the modes and frequencies 
of Im-vacuo ring-stiffened cylindrical shell. Beeanse of the simple 
extension required to apply these frequencies to the quagi-tmpedance 
developed, the general form of an approximate metaod by Ble ich! 45) 
as reported in (43) are included here. Results of the eporoximeie 
engineering analysis of (45) yield resuits in acreement with model 


experiments within + 10 percens. 


Nomenclature 
W. = frequency of shell 
Dy = frequency of ring 
a, = mass of shell per unit of area 


Fy mags of ring frame per unit leagth 
St = frame specing 
x = gecond moment of inertia of ring frame 


| 


section plus a length of shell plating 


equal to 2/p. 


5 


K& K <= Bleich constants tebulated in (45) 
d = thickness of shell 


= Poisson's ratio 


4 
it 


E = Young's Motilus 


a = mean shell radius 
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Bleich represents the total frequency of the in-vacuo 
vibratios of stiffened cylinders es a summation of the shell 


frequencies plus the ring frequencies, 








2 2 2 
O =O + iy (D~-1) 
2 2 2 
0” .* fa) a a x 2B (De2) 
2 m 
a’ +1 a r 
2 -2 2 
ag = Be x —x{i+&x%) (D-3) 
a 2(l-v) m, K. a 
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ni «at xm _— 3 
8 = = 
oy “° *y 
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2 
al [tr+$ (n® - 2) - 2 (a? ~ 1) 6 (x + plan 
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(2) 
(2) 
(3) 
(4) 


(5) 


(6) 
(7) 
(3) 
(9) 
(10) 


(22) 
(12) 


(13) 


(24) 


(15) 
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